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ABSTRACT 

We analyze the vertical distribution of the resolved stellar populations in six low-mass (V max = 
67 — 131 km s _1 ), edge-on, spiral galaxies observed with the Hubble Space Telescope Advanced 
Camera for Surveys. In each galaxy we find evidence for an extraplanar stellar component extending 
up to 15 scale heights (3.5 kpc) above the plane, with a scale height typically twice that of 2-D fits to 
K s band 2MASS images. We analyze the vertical distribution as a function of stellar age by tracking 
changes in the color-magnitude diagram. The young stellar component (<10 8 yrs) is found to have 
a scale height larger than the young component in the Milky Way, suggesting that stars in these low 
mass galaxies form in a thicker disk. We also find that the scale height of a stellar population increases 
with age, with young main sequence stars, intermediate age asymptotic giant branch stars, and old 
red giant branch stars having succesively larger scale heights in each galaxy. This systematic trend 
indicates that disk heating must play some role in producing the extraplanar stars. We constrain 
the rate of disk heating using the observed trend between scale height and stellar age, and find that 
the observed heating rates are dramatically smaller than in the Milky Way. The color distributions 
of the red giant branch stars well above the midplane indicate that the extended stellar components 
we see are moderately metal-poor, with peak metallicities around [Fc/H]= —1 and with little or no 
metallicity gradient with height. The lack of metallicity gradient can be explained if a majority 
of extraplanar RGB stars were formed at early times and are not dominated by a younger heated 
population. Our observations suggest that, like the Milky Way, low-mass disk galaxies also have 
multiple stellar components. In its structure, mean metallicity and old age, the RGB component in 
these galaxies seems analagous to the Milky Way thick disk. However, without additional kinematic 
& abundance measurements, this association is only circumstantial, particularly in light of the clear 
existence of some disk heating at intermediate ages. Finally, we find that the vertical dust distribution 
has a scale height somewhat larger than that of the main sequence stars. 

Subject headings: galaxies: spiral galaxies: structure galaxies formation galax- 

ies:individual(IC 5052, NGC 55, NGC 4144, NGC 4244, NGC 4631, NGC 5023) - 
dust 



1. INTRODUCTION 

Galactic structure has long been recognized as a key 
constraint on theories of galaxy formation. The presence 
of two distinct components - disks and spheroids - sug- 
gests that at least two separate physical mechanisms were 
active during the formation of spiral galaxies. Within the 
Milky Way, however, detailed analyses of the kinematics 
and chemistry of nearby stars have revealed additional 
components: a thin disk with scale heights dependent on 
age, a thick disk and a stellar halo. Each of these compo- 
nents place unique constraints on the Galaxy's evolution. 
In particular, the ages and metallicities of the old thin 
disk, thick disk and halo indicate that they are remnants 
of the initial stages in the assembly of the galaxy and 
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thus constrain the early evolution of galaxies. 

Several scenarios for the creation of the old thin and 
thick disks are consistent with current observations in the 
Milky Way. These scenarios can be categorized into three 
main types: (1) the creation of a thick disk from a thin 
disk by heating by molecular clo uds, spiral arms, star for- 
mation or accretion events (e.g.lSpitzer fc Schwarzschild l 
Il95l } iBarbanis fc Woltrerl 11967] lLacevI 119911 iKrounal 
120021 lOuinn et alJll993HGnedinll2003D . (2) the slow col- 
lapse of the proto- Galaxy forming the thick and thin disk 
in succession (e.g. lEggen et al.l ll962: Gilmore 1984J), and 
(3) the formation of a thick disk from mergers, either by 
direct accretio n of stars or by in situ formation from ac- 
creted gas (e g.lBekki fc Chibal200 lllGilmore et al.l2002t 
lAbadi et all 120031 iBrook et al] 12004ft . In this last sce- 
nario the thin disk would likely form by a settling of 
gas from the merger events. It is widely assumed that 
the stellar halo forms from accreted satellites (see review 
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bv lFreeman &: Bland- Hawthornl2002[) . Determination of 
which processes are important in galaxy formation re- 
quires the study of stellar components in other galaxies. 

Unfortunately, detailed analyses of older stellar com- 
ponents are difficult outside the Milky Way. Some thick 
disks and stellar halos have be en identified using broad- 
band surface photometry (e.g. |Bursteinl 119791: iTsikoudil 
19791 iDalcanton fc Bernstein! 120021 IPohlen et all 120041 
Zibetti et alJl2004|) . but their low surface brightness pre- 
cludes all but a most cursory study of their structure 
and stellar content. Recently, HST imaging has begun 
to allow richer analyses of the resolved stars in the thick 
disk and hal o l|Brown et al.ll2003 iTikhonov et~aTll2005t 
lMouldll2005|) . The resulting color -magnitude diagrams 
can provide robust constraints on the metaUicities and 
ages of the stars which dominate these old components 
in other galaxies. 

In this paper we take advantage of a large database of 
ACS imaging of nearby edge-on galaxies to analyze the 
stellar populations of the galaxies as a function of height 
above the galaxy midplane ("disk height"). The wide 
field of view of ACS allows us to track changes in the 
color-magnitude diagram to more than 3 kpc above the 
plane. Moreover, our sample galaxies have much lower 
masses than the Milky Way, and thus constrain how the 
extraplanar stellar content varies with galaxy mass. The 
properties of the sam ple and thei r globa l color-magnitude 
diagrams are given in S eth et alJ l)2005|) (Paper I). In this 
paper, we study the spatially resolved stellar populations 
in eight fields of six galaxies that were close enough to de- 
termine accurate Tip of the Red Giant Branch (TRGB) 
distances. In particular, we show that there is a signif- 
icant stellar population high above the plane of these 
galaxies and that these stars are old and metal-poor, 
comparable to high-latitude stellar populations in the 
Milky Way. 

In §2, we review the results of Paper I for these six 
galaxies and describe the completeness corrections used 
to correct the stellar number counts in the remainder 
of the paper. In §3 we show that we trace stars in the 
host galaxies out to many scale heights and that these 
stars have a larger scale height than expected from fits 
to 2MASS K s band images. In §4 we select three popu- 
lations of stars from color-magnitude diagrams (CMDs) 
and show that the scale height of a stellar population 
increases with age. We derive metallicity distribution 
functions for the extended stellar component in §5 and 
then conclude by discussing our results in a broader con- 
text in §6. 

2. REVIEW OF THE GALAXY PROPERTIES 

Table 1 shows the position, type, maximum circular 
velocity (V mSL x), distance modulus (m — Mtrgb), scale 
length (h r ) and the K s band half-light height (z 1 / 2 ) of 
the six galaxies we will discuss in this paper. The lat- 
ter two parameters were determined from 2-D model 
fits to 2MASS K s band data, presented in Paper I. 
The vertical component of these models is defined us- 
ing the distribution of an iso thermal population of stars 
(van der Kruit fc SearlellT98l : 

£(z) oc sech 2 (j^J (1) 



where S(z) is the surface brightness or density at a po- 
sition z above the midplane, and zq is the scale height. 
We will use this functional form to fit vertical distribu- 
tions throughout this paper. Note, however, that this is 
one of many equations commonly used to describe the 
vertical distribution of star s in galaxies, a goo d overview 
of which can be found in ijPohlen et alJ l200(l). Most of 
these functional forms vary near the midplane, but have 
similar exponential declines at large disk heights. When 
comparing galaxies in this paper, the disk heights will 
be normalized by the z X j 2 parameter, which gives the 
height containing 50% of the K s band light in the model 
fits. It is related to z , 2^/2=0.549 z , and is similar 
to the exponential scale height h?., which at large scal e 
heights is equal to \zq ijvan der Kruit fc Searld 11981). 
The values of z X j 2 for the six galaxies range from 160- 
280 pc. For comparison, the Milky Way thin disk has 
expon ential scale he ights rangi ng from ^100 pc at young 
ages i|Schmidtll963|) to 330 pc (jChen et al.l200H) at older 
ages. 

All six galaxies are within 8 Mpc and are type Sc or 
later. The maximum circular velocities are all below 135 
km/sec suggesting that these objects are closer in mass 
to the LMC than to the Milky Way. We note that all the 
galaxies except NGC 4631 have circular velocities well 
below 120 km/sec, which appears to mar k a break in 
the dust pr operties llDalcanton et al.ll2004[) and current 
metallicity l|Garnettll2002fl in spiral galaxies. The scale 
lengths of our sample galaxies range from 0.9 to 1.6 kpc, 
^2-3 times smaller than the Milky Way scale length 1 . 
None of the galaxies has an apparent bulge component, 
although NGC 4244 does have a prominent central stellar 
cluster, which is clearly visible in the ACS and 2MASS 
K s band images. 

Our observations include eight HST/ ACS fields in the 
six galaxi es shown in Tab le 1 and are described fully in 
Paper I ijSeth et al.ll2005|) . These observations were ob- 
tained to study the dust-lane properties of these galaxies, 
and hence are centered on the galaxy midplane. The di- 
mensions of these eight fields are given in Table 2, which 
shows the minimum and maximum disk radius in terms 
of the scale length and the minimum and maximum disk 
height in terms of z 1 / 2 . In general, each field is lo- 
cated close to the center of the galaxy. However, two 
of the galaxies, NGC 55 and NGC 4631 have additional 
fields located further out in the disk that are given a 
'-DISK' suffix. Note that many of the galaxies lie diag- 
onally across the chip meaning that the extremities of 
the ranges given in Table 2 are not well sampled. In this 
paper, we focus on the vertical distributions of the stars, 
and where not otherwise noted, analyze all data at the 
same disk height together. 

The approach proposed above is valid as long 
as the scale height of disk components does not 
vary substantially with radius, an assumption that 
has been v erified through observations of edge-on 
galax ies (e.g. Ivan der Kruit fc Searielll981l IPohlen et alJ 
l2004[) . We note however, that there are some anal- 
yses which indicate that the scale height of galaxies 

1 The Milky Way thin & thick disk scale lengths are r ather un- 
certa i n with recent values ranging between 2- 4 kpc (e.g. INg et all 
1199a iMendez k, van Altencjll998l:IOihcJI2001l) . we will use a scale 
length of 3 kpc for comparisons to the Milky Way made in this 
work. 
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TABLE 1 

Galaxy Sample Properties 



Galaxy 


RA(J2000) 


Dec(J2000) 


incl. 


Type 


km/sec 


m - M TRGB 


h R 
["] 


h R 
[kpc] 


z l/2 

H 


z l/2 

[pc] 


IC 5052 


20 52 02.9 


-69 11 45 


89 


SBcd 


79 


28.90 


53.87 


1.57 


7.33 


214 


NGC 55 


00 14 54.4 


-39 11 59 


80 


SBm 


67 


26.63 


93.30 


0.96 


23.34 


240 


NGC 4144 


12 09 58.3 


+46 27 28 


83 


SBc 


67 


29.36 


30.56 


1.10 


6.96 


251 


NGC 4244 


12 17 29.5 


+37 48 28 


84.5 


Sc 


93 


28.20 


84.27 


1.78 


12.15 


257 


NGC 4631 


12 42 07.7 


+32 32 33 


85 


SBcd 


131 


29.42 


35.49 


1.32 


7.54 


280 


NGC 5023 


13 12 11.7 


+44 02 17 


88 


Sc 


77 


29.10 


39.81 


1.28 


5.10 


160 



Note. — RA, Dec, TRGB distance moduli, h r , and Z\i% taken from Paper I. Galaxy Types and V max from 
H YPERLEDA /LED A IPaturel etlT1[l995ll2(}03T). Inclinations tak e n fromlBecker et alj tl98ll : IHummel et alj 
11989) : IMartinl 119981) : IDTTingI 119961) : IHummel fe Dettmarl 119901) : Ide Griis fc van der Krufa 119961) (in same 
order as table). 



TABLE 2 
Field properties 



Field # Stars R min R max z min 2 max 

[R/h r ] [R/h r ] [z/z 1/2 ] [z/z 1/2 ] 



IC 5052 


68636 


-2.9 


1.9 


-22.1 


15.6 


NGC 55 


281536 


-1.6 


1.4 


-5.8 


5.3 


NGC 55-DISK 


253108 


3.3 


6.2 


-7.3 


5.5 


NGC 4144 


60552 


-2.3 


5.9 


-18.7 


14.6 


NGC 4244 


121238 


-1.4 


1.2 


-6.8 


10.5 


NGC 4631 


104940 


-3.1 


3.4 


-17.8 


10.2 


NGC 4631-DISK 


97656 


-10.1 


-1.6 


-21.2 


14.8 


NGC 5023 


42293 


-2.7 


3.1 


-18.1 


25.3 



As described in Paper I, we conducted extensive arti- 
ficial star tests to characterize the completeness in each 
field as a function of magnitude, color and location. The 
galaxies in our sample have high surface brightnesses 
near their midplanc making the completeness there much 
lower (see Fig. 2 from Paper I) . The goal of this paper is 
to analyze the vertical distribution of stellar populations, 
and therefore correcting for this varying completeness is 
critical. 

For each field, over 5 million artificial stars were in- 
serted at random positions, with magnitudes between 
F606W of 18 and 29, and with F606W-F814W colors 
between -1 and 3. These stars were then run through 
the same pipeline used to determine the stellar photom- 
etry. Artificial stars that coincided with actual sources 
were considered detected only if the input magnitude of 
the artificial star was brighter than the actual source in 
both bands. 

To enable completeness corrections for individual stars, 
we determined our completeness for the artificial stars in 
bins of magnitude, color and local surface brightness. For 
the magnitude and color, we used 0.f5 magnitude wide 
bins. At its steepest, the completeness as a function of 
magnitude varies by at most 6% over 0.15 magnitudes, 
so any error introduced by the binning should be smaller 
than that. We then determined the size of the surface 
brightness bins such that there would be ~50 stars in 
each final bin. Determining the completeness from 50 
stars gives a random error in the completeness of ~6%. 
An aperture around each star from 11 to 14 pixels was 
used to determine the local surface brightness. We deter- 
mined the surface brightness level from the F606W image 
so as to include the effects of the HII region emission visi- 
ble in a number of galaxies (most notably NGC 55, which 
therefore has the brightest completeness level in Fig. 2). 
This emission is not seen at all in the F814W images, due 
to the lack of strong emission lines in F814W bandpass. 
Using this binned completeness function, we are able to 
determine the completeness of any individual star based 
on its color, magnitude and local surface brightness level 
to within ~10%. In the stellar density profiles presented 
in §3 & §4, the completeness corrections are up to 200% 
near the midplane, but fall to < 30% at z/ zi/ 2 > 3. 

In addition to correcting for the completeness, magni- 
tude limits must also be set to insure that we are not us- 
ing stars fainter than we can detect in the higher surface 



flares with increasing radius, both in our own Galaxy 
I Lopez-Corredoira et al.l l2002fl and in edge-on galaxies 
I de Griis fc Peletieilll997t iNaravan fc ,Iodl2002l) . How- 
ever. Ide Griis fc Peletierl l)1997() show that in their sam- 
ple of edge-on galaxies, late-type galaxies such as those 
observed here have little or no flaring. In one of our 
galaxies, NGC 4244, flaring of the H I gas by a fa ctor of 
^3 is seen between radii of 8-13 kpc (Oiling 1996). Even 
if similar flaring occurs in the stellar distribution, the 
star counts at the disk heights we consider here will still 
be dominated by stars located (radially) near the center 
of the galaxy. We note that we do see some evidence for 
modest flaring in our two '-DISK' fields which lie at large 
radii (see §4.2.1). 

In Paper I we presented the color-magnitude dia- 
grams (CMDs) for each of the galaxy fields discussed 
here. These are reproduced in Figure 1 which shows 
the F606W-F814W color vs. the F814W magnitude. 
All magnitudes in this paper are given in the VEGA- 
mag system (see Paper I for details) . Each of the CMDs 
show plumes at F606W-F814W colors of -0.1 and -1.3 
which are young main sequence (MS) and helium burn- 
ing (HeB) stars respectively. An old RGB can be seen 
in each galaxy at colors of —1 and an intermediate-age 
asymptotic giant branch (AGB) extends from the tip of 
the RGB to redder colors. Further details on the com- 
ponents can be found in §4.1 and in Paper I. The boxes 
shown in Figure 1 isolating these components will be dis- 
cussed in greater detail in §4. 

2.1. Completeness Corrections 
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Fig. 1. — Color-magnitude diagrams of each galaxy. Contours are used where the density of points becomes high. The contours are 
drawn at densities of 75, 100, 150, 200, 250, 350, 500, 750, 1000, and 1500 stars per 0.1 magnitude and color bin. The bottom-left panel 
contains a key to the gray lines which demarcate sections of the CMD. The line across the bottom shows the 20% completeness limit in 
bright regions of the galaxy (see §2.1 for details) and the other lines show the regions used in each galaxy for the MS/HeB, AGB and RGB 
stars. 



brightness regions of the image (i.e. the midplane). We 
therefore choose to limit our analysis to regions of the 
CMD that fall above a conservative 20% completeness 
limit in regions of high surface brightness. Figure 2 shows 
the 20% completeness limits for the brightest regions in 
each field. As can be seen, the completeness limit rises 
towards redder colors, and steepens at colors redder than 
F606W-F814W of 1. To make a smooth boundary that 
is easily applied to our data, we fit the 20% completeness 
curves to two lines intersecting at F606W-F814W=1. Ta- 
ble 2 shows the results of these fits, by giving the com- 
pleteness limit at F606W-F814W of -1 (F814Wi im _i), 1 
(F814Wi imj i), and 3 (F814Wi imj3 ). We use these'limits 
throughout this paper to insure that comparisons made 
between stellar populations at different disk heights are 
valid. 

We note here that although we can correct for incom- 



pleteness due to crowding, we cannot correct for the at- 
tenuation of stars by dust. We will show in §4 that all 
the galaxies in our sample are optically thick near their 
midplanes. Therefore at low galactic latitudes, our com- 
pleteness corrected stellar census will fall short of the 
true number of stars. 

3. VERTICAL DISTRIBUTION OF STARS 

We demonstrate in this section that there are signif- 
icant numbers of stars well above the planes of all our 
disks and that the profiles of these stars do not fit the 
profiles expected from the ground-based K s band galaxy 
fits. 

In Figure 3 we present the completeness corrected sur- 
face density profiles of all the detected stars (£ a n) above 
the completeness limits given in Table 3. Two lines are 
shown for each galaxy, giving the profile on both sides 
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Fig. 2. — The 20% completeness limit in the brightest sections 
of each field as a function of color. The variation from ficld-to-field 
is due to increased crowding / higher surface brightness regions in 
some fields. 



TABLE 3 
20% Completeness Limits 



Field 


F814W lim ,_! 


F814Wi imil 


F814W lim , 3 




[mag.] 


[mag.] 


[mag.] 


IC 5052 


26.48 


26.05 


24.51 


NGC 55 


25.39 


24.86 


23.51 


NGC 55-DISK 


26.52 


25.89 


24.28 


NGC 4144 


26.72 


26.28 


24.67 


NGC 4244 


26.04 


25.65 


24.07 


NGC 4631 


25.94 


25.49 


23.88 


NGC 4631-DISK 


26.30 


25.83 


24.18 


NGC 5023 


26.76 


26.52 


24.88 



10.00 



-i — i 1 — i 1 — i r- 

1C5052 
NGC0055 
NGC0OS5-D1SK 
NGC4144 
NGC4244 
NGC4631 
NGC4631-DISK 
NGC5023 




10 

\x\M 



15 



20 



Fig. 3. — The surface density of stars as a function of height 
from the midplane. Curves have been completeness corrected as 
described in the text. Each galaxy has two lines, one for above 
and below the plane. The dashed lines indicate stellar distributions 
with scale heights 1 & 2 X that of the K B band fits to the disks. The 
dot-dashed lines for the NGC 4631 fields indicate the regions where 
contamination of stars from companion NGC 4627 is significant. 



of the disk. To determine the surface density, stars were 
binned as a function of scale height. After binning, the 
completeness-corrected number of stars was divided by 
the area of the bin to obtain the surface density. Only 
bins with an area of more than 300 arcsec 2 are plotted. 

Figure 3 shows that we trace the stellar component 
of each galaxy out to large disk heights, with several 
galaxies being traced out beyond 10zi/ 2 . The profiles 
in general are fairly symmetric, the most notable excep- 
tion being NGC 4631 and NGC 4631-DISK. These fields 
are contaminated on one side (shown using dot-dashed 
lines) by the presence of the companion galaxy NGC 4627 
(see Fig. 3 from Paper I). The decrease of the profiles 
with increasing scale height out to the edge of the fields 
strongly suggests that the profiles remain above the sur- 
face density of foreground Galactic stars and background 
unresolved dwarf galaxies in our magnitude range. Only 
above 10zi/2 in IC 5052 and NGC 5023 is there some 
evidence for the leveling off that would be expected as 
we reach the foreground/background level. Note that the 
galaxies are all at galactic latitudes above 35°. We there- 
fore can safely assume that a vast majority of stars in our 
images are located in the host galaxies and we make no 
corrections for foreground/background sources. 

The dashed lines in Figure 3 show sech 2 profiles with 
scale heights one and two times the measured K s band 
scale height (note that because the plot is scaled by the 
K s band scale height these profiles are the same for each 
galaxy). As described in detail in Paper I, models were 
fit using a Levenberg-Marquardt algorithm with uniform 
weighting on all unmasked pixels to K s band 2MASS 
data. Because the 2MASS surface photometry has a lim- 
iting isophote of K s ~ 20.0 mag/arcsec 2 , only relatively 
high surface brightness features could be fit. The typical 
K s band peak surface brightness of our galaxies is ~18 
(Paper I, Table 2), which means that galaxies are only 
detected in the K s images out to a few zy2 from the 
midplane. Although the K s band light is often thought 
to be dominated by the RGB stars that trace an old stel- 
lar population, we will show in §4 that in these low mass, 
late-type galaxies, it more closely traces the young and 
intermediate age populations and is thus dominated by 
red supergiant and AGB stars. 

Figure 3 shows that the outer portions of the stellar 
density profiles of the galaxies appear to be broader than 
indicated by the K s band scale height. To quantify this, 
we fit the stellar density profiles between 5ziy 2 and 10zi/2 
on each side of the midplane in each galaxy. Only profiles 
with data beyond 8zi/ 2 were fit (thus excluding NGC 55, 
NGC 55-DISK and one side of NGC 4244). In the Milky 
Way, the stellar profile deviates f rom the thin disk profil e 
beyond ~1 kpc above the plane ijGilmore k. Reid| [1983). 
For our galaxies, 5zi/ 2 is ~1 kpc. Therefore we would 
expect our fitting range to be sensitive to a possible thick 
disk component in these galaxies. The scale heights of 
these fits above and below the midplane are shown in 
the third and fourth columns (respectively) of Table 4. 
All of the fitted stellar profiles are significantly broader 
(x 1.5-2.4) than the K s band scale height. This observa- 
tion strongly suggests that these galaxies contain a more 
broadly distributed stellar component not traced by the 
K s band 2MASS images. 

We note that the fitted scale heights, both in the K s 
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TABLE 4 

Resolved Stellar Component Scale Heights 



Field 


zq [pc] 


All + 
z Q [pc] 


All - 

2 [pc] 


MS 
zq [pc] 


AGB 
ZQ [pc] 


RGB 

z [pc] 


fir / Zn I\/TQ 


fir / Zn APR 


fir I Zi\ RC R 


IC 5052 a 


390 


767±30 


686±44 


261±11 


467±6 


655±6 


6.0 


3.4 


2.4 


NGC 55 


437 






327±24 


644±10 


701±3 


2.9 


1.5 


1.4 


NGC 55-DISK 


437 






526±1 


741±15 


999±6 


1.8 


1.3 


1.0 


NGC 4144 


457 


940±41 


965±15 


374±27 


699±16 


934±18 


2.9 


1.6 


1.2 


NGC 4244 


468 


740±40 




325±20 


443±24 


551±9 


5.5 


4.0 


3.2 


NGC 4631 b 


510 


927±46 




510±26 


895±51 


1154±194 


2.6 


1.5 


1.1 


NGC 4631-DISK b 


510 


1131±50 




505±22 


1200±1 


1387±73 


2.6 


1.5 


1.0 


NGC 5023 a 


291 


505±32 


534±26 


204±6 


289±6 


391±4 


6.3 


4.1 


3.3 



Note. — (a) In IC 5052 and NGC 5023, MS, AGB and RGB fits excluded disk heights below 1.5z 1/2 . All other galaxies 
excluded below 3z 1 /2- (b) Negative values of z excluded beyond 2 kpc due to the presence of companion galaxy NGC 4627. 



band and using stellar density profiles, can differ from the 
true scale height due to a number of factors. First, dust 
attenuation can obscure light near the midplane of the 
galaxy. In the K s band this attenuation should be small, 
and we ameliorate this problem in the fits presented here 
by avoiding the midplane. Second, the galaxies may not 
be perfectly edge-on. Based on previous observations, 
the least inclined of the galaxies is NGC 55, which has 
a -80° inclination l|Hummel et alJ ll98fi1. This would 
give a fitted scale height ^30% greater than the intrin- 
sic disk scale height. We also note that NGC 55 and 
NGC 4631 are fairly irregularly shaped making the fits to 
these galaxies less reliable than for the other four galax- 
ies. 

4. VARIATION IN DISTRIBUTION WITH STELLAR 
POPULATION 

We now turn our analysis to stars selected in regions of 
our CMDs that isolate stellar populations with different 
ages. Using this method we show that the older stellar 
populations have an increasing scale height. In §4.3, we 
examine the variation of scale height in the context of 
disk heating models. We then present simplistic dust 
models in §4.4 and show that the dust extinction in these 
galaxies is distributed in a component that is broader 
than the young stellar populations. 

4.1. Selection of CMD regions 

We attempt to separate our data into young, interme- 
diate, and old stellar populations by selecting stars from 
different regions in the CMD. The young stars are found 
in the Main Sequence (MS) component and in the red 
and blue Helium Burning (RHcB,BHcB) sequences (see 
Paper I, Fig. 1 for a schematic CMD), all of which should 
contain stars under a few 100 Myr in age. For the inter- 
mediate age stars we select AGB stars brighter than and 
redward of the RGB, resulting in ages ranging between 
a few 100 Myr and a few Gyr. Lastly, for the old pop- 
ulation of stars we select RGB stars, which have ages in 
excess of 1 Gyr, although some AGB stars will also be 
found in the same region. 

The actual regions used for the selection are shown in 
the CMDs in Fig. 1. The bottom right figure is a cartoon 
illustrating the selection of the MS/HeB, AGB and RGB 
regions. The RGB region was selected using lines with 
slopes of 3.3 and 6.6 and F606W-F814W colors of 1.0 and 



1.6 at the tip of the RGB. The MS region was defined 
by taking all stars redward of -0.5 and blueward of a 
line with slope of 3.3 and a color of 0.7 at the TRGB 
magnitude. Finally, the AGB region isolates stars less 
than two magnitudes brighter than the TRGB magnitude 
and redwards of a line with slope of 3.3 and a color of 
1.2 at the TRGB magnitude. These boundaries were 
combined with the TRGB magnitude given in Table 1 
and the completeness limits from Table 3 to determine 
the final regions for each galaxy shown in Figure 1. The 
regions were chosen somewhat conservatively - e.g. we 
chose to put space between the MS /HeB section and the 
RGB so that there would be little overlap between the 
two due to dust extinction or large photometric errors at 
faint magnitudes. From here we will refer to the stars 
in these CMD boxes as the MS, AGB and RGB stellar 
populations. 

4.1.1. Synthetic CMDs 

To determine the typical ages of stars detected in 
our CMD boxes and to facilitate quantitative com- 
parisons between galaxies and their different stellar 
populations, we ge nerated synthe tic CMDs using the 
MATCH program iDolphinl l20f)l) and the IAC-STAR 
program 2 (Aparicio & Gallart 2004), using isochrones of 
iBertelli et all l)1994|) and iGirardi et alJ l|2000|) in both 
cases. The synthetic stars were generated assuming a 
constant star formation rate (SFR) from 13 Gyr ago 
to the present, and a metalli city that steadi ly increased 
from [Fe/H]= -1.7 to -0.4 ilCnrnettlEool. We used 
slightly different IMFs in the two CMDs. For the 
MATCH CMD a pure Salpeter IMF {a = 2.35) is as- 
sumed between 0.1 and 120 M , whereas in th e IAC- 
STAR CMD we used the default iKroupa et all (1993) 
IMF, which is steeper at the high mass end (a — 2.70). 

To compare these CMDs to our observations, the syn- 
thetic stars were first transformed from Johnson V & I to 
VEGAmag F606W & F814W colors. We then mimicked 
observations of each galaxy as follows. First, each star 
was randomly assigned a surface brightness value based 
on the values of detected stars in each frame. Then, 
using the artificial star tests, we determined the chance 

2 This work has made use of the IAC-STAR Synthetic CMD 
computation code, IAC-STAR is supported and maintained by the 
computer divison of the I AC. 
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each star was detected and a magnitude error based on 
the star's initial F606W & F814W magnitudes (assuming 
the distance moduli shown in Table 1) and the surface 
brightness value. A final CMD was then made by ran- 
domly determining if each star was detected and applying 
the determined errors. The resulting CMDs looked qual- 
itatively similar to our observed CMDs, with the most 
notable difference being an offset of the AGB stars to 
somewhat brighter magnitudes in both synthetic CMDs 
relative to the real data and a deficit of MS/HeB stars 
in the 1AC-STAR CMDs relative to the MATCH and 
the real galaxy CMDs. This could indicate either that 
the galaxies have enhanced rec ent star formation or that 
their IMF is not as steep as the lKroupa et alJ l)1993|) IMF 
on the high mass end. 

Figure 4 shows the resulting age distribution of the 
MS, AGB, and RGB boxes in NGC 4144 using the 
MATCH and IAC-STAR synthetic CMDs. The age dis- 
tributions for other galaxies are similar. This figure 
clearly demonstrates that we are separating the stars 
into young, intermediate-age and old populations with 
our chosen CMD boxes. However, the separation is not 
perfect. Each bin has significant overlaps with the others 
due to unavoidable photometric errors and to true over- 
lap in the colors and magnitudes of stellar populations 
with different ages. The IAC-STAR CMDs have age dis- 
tributions similar to the MATCH CMDs, but with the 
AGB populations weighted more towards older ages and 
a more significant contamination of old stars in the MS 
box (probably due to the relative lack of MS stars in the 
IAC-STAR models). Both effects likely result from the 
steeper IMF assumed for the Aparicio CMD. In the fol- 
lowing sections of the paper, we use the MATCH CMD 
for comparisons with observations because it more closely 
reproduces the ratio of young MS and HeB stars relative 
to the number of older stars. 

We note that these sunthetic CMDs assume a constant 
star formation rate and thus are not useful in determining 
true star formation histories. However, we will be able to 
use them to get a sense of relative star formation histories 
(SFHs) as a function of scale height and to get a rough 
sense of the ages of the stars in our CMD regions. 

4.2. Stellar Density Profiles 

Now we compare the surface density profiles of the 
MS, AGB and RGB stars to examine possible varia- 
tions in stellar population with disk height. Figure 5 
shows the completeness-corrected profiles as a function 
of disk height for each field. Each profile is derived using 
the same methodology as in §3, typically using ^10,000 
stars per field. The surface densities are then normalized 
to have J J^dz = 1. All the fields show a similar pat- 
tern. The MS (blue) stars have the narrowest distribu- 
tion while the AGB (red) and RGB (orange/yellow) stars 
have broader distributions and typically show a dip near 
the midplane. Because we have corrected for incomplete- 
ness, the dip almost certainly due to dust absorption, as 
we demonstrate in §4.4 with a very simple model. 

Figure 5 suggests that older stellar populations become 
more prominent with increasing disk height. We quan- 
tify this trend in Figure 6, which shows the ratios of 
surface densities in our different age bins. The ratios 
were normalized to those expected for a constant SFR 
using the MATCH synthetic CMDs (see §4.1.1). A ratio 



of one in Fig. 6 therefore corresponds to a constant star 
formation rate and increasing values correspond to older 
stellar populations. We note that the ratio is only plot- 
ted where the signal-to-noise of the ratio is greater than 
3. The small number of MS stars at large scale heights 
limits our ability to trace the Ergb /Ems and Eagb/Ems 
as high above the midplane as the profiles shown in Fig- 
ures 3, 5 and 7. Also, NGC 4631 and NGC 4631-DISK 
are not included in Figure 6 because the high complete- 
ness limit results in very few RGB stars (see Fig. 1) and 
an increased contamination of AGB stars in the RGB 
box. 

The top and middle panels of Figure 6 show that in 
each of the fields, RGB stars become more numerous rel- 
ative to MS and AGB stars with increasing disk height. 
However, this trend shows an enormous variation from 
galaxy to galaxy. In IC 5052 the ratio Ergb/Ems be- 
comes as high as ~100 times the midplane value, while 
in NGC 55 and NGC 4244 the increase is much more 
moderate, to <10 times the midplane value. This vari- 
ation is most likely the result of a range of recent SFRs 
in our galaxies. The ratio Ergb/Eagb is much more 
consistent from galaxy-to-galaxy, however. This may re- 
sult from the overlapping time range spanned by stars in 
the AGB and RGB boxes and/or the large time ranges 
these boxes span relative to the MS. We argue in §5 that 
the RGB population is likely to be dominated by truly 
old stars much older than the AGB population. Inter- 
estingly, the field showing the the flattest Ergb/Eagb 
profile was the NGC 55-DISK field located in the outer 
parts of the NGC 55 disk, perhaps suggesting a differ- 
ent star formation or dynamical history at large radii 
(~5x ha). However, this galaxy is the least inclined in 
our sample and is somewhat irregular in shape, therefore 
results for this one system should not be overinterpreted. 

The low values of Eagb/Ems and the high values of 
Ergb/Eagb result from a lack of AGB stars compared to 
the constant SFR MATCH synthetic CMD. This would 
seem to suggest that the galaxies' star formation histo- 
ries (SFHs) are depressed at intermediate ages and en- 
hanced at young ages. However, as we noted in §4.1, the 
AGB morphologies in the synthetic CMDs are not well 
matched to the observational CMDs, probably because 
of the difficult y in modeling the AGB phase of evolution 
<IMarigoll200ll) . This discrepancy combined with the dif- 
ferences seen between the two sets of synthetic CMDs 
suggests that a derivation of accurate SFHs using just 
the brightest stars in a galaxy is not yet possible. 

To check if the scatter in the Ergb/Ems ratio was in 
part due to varying radial coverages of the galaxies (Ta- 
ble 3), we remade the plots in Figure 6 using only stars 
within the central scale length of each galaxy. These 
plots were similar to those shown and showed compara- 
ble scatter. This suggests that the observed variations 
from galaxy to galaxy in the stellar populations ratio 
reflect global differences in the galaxies' SFHs and/or 
vertical structure. For instance, if we assume the trend 
towards older populations with increasing scale height re- 
sults from disk heating, then the scatter in Figure 6 sug- 
gests substantial variations between galaxies in either the 
mechanisms that heated the disk, or the SFH of the disk. 
Despite these variations, Figure 6 gives strong evidence 
that overall the age of the stellar populations increases 
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Fig. 4. — Histogram of ages detected in denned CMD boxes (see Fig. 1) for NGC 4144 assuming a constant star formation rate from 
13 Gyr to the present. Histograms are based on synthesized CMDs created with the MATCH (left) and IAC-STAR (right) programs as 
described in §4.1. These plot shows that stars in the MS box (blue) are dominated by stars ~100 Myr in age, while the AGB (red) and 
RGB (orange/yellow) boxes have typical ages of ~1 Gyr and ~10 Gyr respectively. Similar plots for other galaxies are qualitatively similar. 
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Fig. 5. — The normalized surface density as a function of scalchcight for young MS (blue), intermediate-age AGB (red) and old RGB 
(orange/yellow) stars. Each surface density distribution was normalized to integrate to one. Note that in all cases the MS distribution is 
the most peaked while the RGB distribution is the widest. 
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Fig. 6. — The ratio of the surface density of different stellar pop- 
ulations as a function of disk height: RGB/MS (top), RGB/AGB 
(middle), and AGB/MS (bottom). Points were normalized based 
on the number counts of a constant star formation rate model, e.g. 
a value of 1 suggests a constant star formation rate, while values 
>1 suggest a decrease in the star formation rate with time. Note 
that the MATCH synthetic CMDs were used in this normalization 
(§4.1.1). Only points with a signal-to-noise > 3 are plotted. The 
plots clearly show that the stellar populations become older with 
increasing disk height. 



with increasing scale height. 

4.2.1. Stellar Population Scale Heights 

To further quantify the differences in the vertical dis- 
tribution of the three CMD regions, we fit each surface 
density profile to a sech 2 function in which the normaliza- 
tion, central position, scale height (zo), and background 
level were all allowed to vary. We fit each profile only 
at disk heights >3ziy 2 to avoid the dips near the mid- 
plane, except in IC 5052 and NGC 5023 where we used 
disk heights >1.5zi/2 to allow fitting of the very narrow 
Sms profiles. Figure 7 shows the resulting fits to the 
Ems, Sagb, and Ergb profile of each galaxy, in the top, 
middle and bottom panels respectively. The observed 
profiles are shown as a solid line, while the best fitting 
sech 2 function is shown as a dashed line. The dotted lines 
at ±3z 1/2 (±1.5z 1/2 in IC 5052 and NGC 5023) delineate 
the region excluded from the fit. The error bars on the 
data points are used to weight the sech 2 fits and reflect 
Poisson errors in the number counts, but do not include 
uncertainties in the completeness corrections. The scale 
height of the best-fitting sech 2 function is shown in the 
upper-left corner of each panel, and the error shown is 
scaled by the square root of the reduced \ 2 of the fit. The 



reduced \ 2 values for a majority of the fits were between 
0.8 and 1.3, but were larger for NGC 55 and NGC 4631 
due to their irregular structure. Scale heights and errors 
for all the fits are shown in Table 4. 

We find that in each galaxy, the MS scale height value 
is the narrowest followed by the AGB and then the RGB. 
In all cases the RGB population is significantly broader 
than the AGB, MS and K s band zq values. This re- 
sult strongly suggests the presence of an older component 
with larger scale height. An analysis of the variations in 
scale height of the MS, AGB and RGB populations with 
scale length in each of the galaxies turned up no obvious 
trends. We also identified no trends with galaxy rota- 
tion speed, due to the small range of masses spanned by 
our sample galaxies. We note that in some cases, the z§ 
derived for all stars (§3) is somewhat larger than the zo 
derived for just the RGB stars. This results from the 
lower scale heights used in the fits to the different stellar 
populations - if the lower limit for the fit to all stars is 
reduced from 5zi/ 2 to 3.21/2, the derived zq is less than 
or equal to the RGB star z$ in each galaxy, as expected 
if the total stellar density is well characterized by a com- 
bination of the MS, AGB and RGB components. The 
fits presented here for these different stellar population 
components are in general not extremely sensitive to the 
range of z values used. Varying the lower limit of the fit 
between \.h-bz\i-i typically changed the AGB and RGB 
zq values by less than 10%. The compact MS components 
were more dramatically affected because of the smaller 
number of stars at large scale heights. 

From this analysis it appears that NGC 4631 has the 
"thickest" old component with an RGB scale height of 
~1250 pc, roughly 2.5 times larger than the MS and K s 
band fits. We note that the fits for NGC 4631 were trun- 
cated at large negative disk heights to prevent the con- 
tamination of stars from companion galaxy NGC 4627. 
IC 5052 and NGC 4144 have similar ratios (-2.5) of RGB 
to MS scale heights, while NGC 4244 has the smallest ra- 
tio, with an RGB scale height only 1.7 times that of the 
MS stars. 

There is also evidence for a modest flaring of the stellar 
components between the central and '-DISK' pointings 
of NGC 55 and NGC 4631. The '-DISK' fields are cen- 
tered 4.8 and 6.1 scale lengths (see Table 2) from the 
center of the NGC 55 and NGC 4631 respectively. An 
increase in zq values by a factor of 1.1 to 1.6 is seen for 
all three components in NGC 55 and for the AGB and 
RGB component in NGC 4631. 

The profiles deviate from the fitted sech 2 profile sig- 
nificantly near the midplane. This deviation is almost 
certainly due to dust, we model this effect in §4.4. At 
larger disk heights (>2-3 kpc) there is also a slight over- 
density in the RGB components of IC 5052, NGC 4144, 
NGC 4244 and NGC 5023. These overdensities hint at 
the possible presence of an even more broadly distributed 
old component. In three of these galaxies, the RGB 
fits had elevated x 2 values relative to the MS and AGB 
fits. We estimate that these overdensities occur at a sur- 
face brightness ^F606W ^ 28 mag arcsec -2 assuming a 
luminosity function similar to galactic globu lar clusters 
(jBuonanno et al.|ll994t iKravtsov et alJll99"7^ . However, 
without better knowledge of the background level, it is 
not possible to verify the existence of this component. 
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The scale heights measured in the K s band (Table 4) 
are closest to those measured for the MS and AGB com- 
ponents. Half of the galaxies have K s band scale heights 
closer to the MS value and the other half closer to the 
AGB value. This suggests that the K s band light in 
these galaxies is dominated by relatively young stellar 
populations, probably red sup ergiants and AGB stars (in 
agreement with the findings of lAoki et alJl99l|) . This re- 
sult runs contrary to the common assumption that the 
NIR light primarily traces older stellar populations (e.g. 
iFlorido et aDI^OOlj) . and is significant in that NIR lu- 
minosity is often used as a proxy for stellar mass when 
comparing galaxies of different types and masses. How- 
ever, we note that our K s band scale heights are biased 
towards higher surface brightness populations due to the 
bright limiting isophote of the 2MASS data from which 
they are derived. 

4.2.2. Comparison to Previous Observations 

The results above indicate that there is a systematic 
increase in the vertical scale heights of older stellar popu- 
lations in our sample of low mass, late-type disks. Before 
investigat ing possible origins for these structural differ- 
ences in S14.3I & Sj3 we now compare our measurements 
of scale heights to previous observations of the vertical 
structure of disks. 

The most detailed constraints on the scale heights of 
different stellar populations come from the solar circle 
of the Milky Way. Studies have revealed a complicated 
disk structure, with a young and old thin disk embedded 
within a more extended thick disk. The young thin disk 
is the narrowest of the three, having a scale height of 
z ~ 200 pc, as trace d by stars wit h bright absolute 
magnitudes (My < 3) l)Schmiddll963D . In contrast, the 
scale heights of the young main sequence stars in our 
sample are almost all significantly larger than the Milky 
Way value suggesting that the low mass galaxies in our 
sample form sta rs in a thicker lay e r than the Milky Way, 
consistent with iDalcanton et alJ l|2004l) . The resulting 
axial ratios for our samples' young star forming disks are 
also much thicker as well, with zo/h r — 1.8 — 6.3 (see 
Table 4) for our sample galaxies, versus Zo/h r ~ 15 for 
the young thin disk of the Milky Way. 

The division of the Milky Way's older stellar popula- 
tions into a thin and t hick disk was first introduced by 
iGilmore fc Reidl l)1983i) to explain a break in the number 
counts of F & G stars at ~ 1 kpc. While the need for two 
old disk components was long debated, recent measure- 
ments of systematic q-element enhancement in thick disk 
stars (most recently | G ratton et alJ 120031: jFeltz ing et all 
12003 iMishenina et all2004UBensbv et alJl2005l) strongly 
suggest that the thick disk is indeed distinct from the 

old thin dis k. Recent observations ijChen et alJ l200l|: 

ISiegel et all 12002) give a scale height for the old thin 
disk of zn ~ 600 pc, similar to found in the original 
IGilmore fc Reidl ©83) study. These same studies sug- 
gest that the exponential scale height of the thick disk 
is h z ~ 700 pc (corresponding to z ~ 1400 pc, thinner 
than originally claimed). The Milky Way thick disk is 
therefore roughly twice the height of the old thin disk, 
and 7 times the height of the young thin disk. 

Within our own sample, the scale height of the old 
RGB component is mostly intermediate between the 
Milky Way old thin disk and the thick disk. Our sample 



galaxies have much lower masses and surface densities 
than the Milky Way, and, lacking any firm model that 
predicts how the properties of the old thin disk and the 
thick disk should vary with galaxy mass, we are hesis- 
tant to attribute the extraplanar population to either an 
old thin or a thick disk on the basis of the the surface 
brightness profiles alone. There are no dramatic inflec- 
tion points in the RGB surface density profiles plotted in 
Figure 7 that would assist in a unique separation of old 
thin disk and thick disk stars, and the possible overden- 
sity of stars above 2-3 kpc may well be due to a stellar 
halo. Even if the RGB component is similar to the Milky 
Way thick disk, this lack of inflection is not unexpected. 
In the Milky Way, the inflection point in the surface den- 
sity of F & G dwarves that marks the separation of the 
thin and thick disks is likely the result of two different 
populations of stars separated in age. The lack of sim- 
ilar inflection points in our RGB profiles can easily be 
explained if the RGB stars don't have as wide a range of 
ages as the Milky Way dwarves. We will show in §5 that 
the RGB stars in our galaxies may very well be domi- 
nated by a single-age population. 

We do note, however, that the axial ratios of the RGB 
disks range between h r / zq — 1.0 — 3.3, with a median of 
1.8 (adopting the ii's-band radial scale length, and av- 
eraging the two independent measurements for NGC 55 
and NGC 4631). For comparison, the axial ratios of the 
old thin and the thick disks of the Milky Way are 5.0 and 
2.1, respectively (assuming h r = 3 kpc for both compo- 
nents). Thus, in terms of their overall structure, the 
RGB component we detect is significantly more round 
than the Milky Way's old thin disk, and is distributed 
more like the Milky Way thick disk. However, without 
additional information we cannot ascribe a common for- 
mation scenario to our observed RGB component and 
the Milky Way thick disk. We revisit this issue in the 
discussion (50, after analyzing the disk heating and the 
vertical metallicity gradients of our sample galaxies. 

Outside of the Milky Way, the most detailed infor- 
mation comes from studies of the vertical distribution 
of resolved stars in HST i mages, similar to the work 
we present in this paper. ITikhonov et alJ l)2005|) and 
ITikhonov fc Galazutdinoval l)2005l) present evidence for 
extended components in six galaxies, which they quali- 
tatively argue correspond to thick disks and halos. Of the 
galaxies that overlap our sample (NGC 55, NGC 4144, 
and NGC 4244), they include archival WFPC2 obser- 
vations to reach greater disk heights in NGC 4244 and 
NGC 55 than spanned by our ACS images. In both 
cases they assume a priori that the RGB stars at 
lower disk heights trace a thick disk. For NGC 4244, 
ITikhonov Sz Galazutdinoval 1)20051) show an exponential 
distribution of RGB stars between ~1 and 3 kpc (their 
Figure 8) that appears to roughly match the scale height 
of the profile shown in Fig. 7. Beyond 3 kpc, they see 
a flattening in the number counts which they claim is a 
halo, but which may also be th e background level. For 
NGC 55. ITikhonov et all l)200 5) plot an exponential dis- 
tribution of RGB stars between 2 and 7 kpc - i.e. at much 
greater disk heights than probed by our data. However, 
based on inspection of their Figure 12, the extended RGB 
component has a zq value of ~ 2 kpc, which is 2-3 times 
the width of the RGB component we fit. Although they 
assume this component is due to a thick disk based on its 



Vertical Distributions of Resolved Stars 



11 




-1500 -1000 -500 500 1000 -4000 -2000 2000 

z [pc] z [pc] 



Fig. 7. — Model fits to the observed surface density distribution as a function of disk height in each galaxy for the MS (top), AGB 
(middle), and RGB (bottom). The dashed line shows the best fitting sech 2 model excluding data with |z/zi/2| < 3 (as shown with the 
dotted lines) while the vertical dotted lines show the range of z values excluded from the fit. The number in the upper left corner gives 
the scale height of the best fitting sech 2 function. 



exponential surface density distribution, the axial ratio 
of this component would in fact be prolate (h r /z Q ~0.5) 
and thus may be more analogous to the Milky Way's stel- 
lar halo 3 . The change in slope also implies a break in the 
RGB distribution in NGC 55 at around 2 kpc. By anal- 
ogy, this may indicate that the marginal overdensities we 
are seeing at comparably large disk heights in our RGB 
profiles might be the signature of an additional broader 
halo component. Mould (200.3) a l so finds the presence of 
old stars at large scale heights, and while these stars are 
automatically assumed to be a thick disk component, no 
detailed analysis of their spatial distribution is presented. 

In addition to these recent studies of resolved stars, 
most previous studies of the vertical structure of 

3 Note, however, that NGC 55 is the least inclined galaxy in our 
sample, complicating the interpretation of its projected structural 
parameters. 



disks have focused on detecting thick disks and stel- 
lar halos using unr esolved surface br i ghtness profiles 
of the galaxies (e.g. iPohlen et all 12004 iFrv et a,lJfl99Ht 
iDalcanton fc Bernstein! 120021: iNeeser et all \200% . Be- 
cause we only detect stars at bright magnitudes, it is 
difficult to accurately convert our measured surface den- 
sity of stars (Figure 2) to a surface brightness. How- 
ever, assuming the outer parts of our galaxies have lu- 
minosity functions sim i lar to Galactic globul ar clusters 
l|Buonanno et all 119941 iKravtsov et all I1997D . we esti- 
mate that we reach F606W surface brightnesses of ^28 
mag arcsec~ 2 . This is comparable to the depth reached 
in deep ground-based observations. 

Only one of our galaxies has been analyzed for ver- 
tically ex t ended components using ground based data. 
IFrv et ail l)1999ft present i?-band surface photometry of 
NGC 4244 and find no evidence of a second thick disk 
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component, based on the lack of an inflection point in the 
surface brightness distribution above the plane. They 
trace the vertical profile of NGC 4244 along the minor 
axis to ~ 2 kpc at which point it falls below their sur- 
face brightness limit of 27.5 mag arcsec -2 . We trace 
the RGB component out to nearly 3 kpc, and find a 
scale height that is similar to their fitted i?-band scale 
height (assuming h z — ^zo). Their lack of an inflec- 
tion hi_the_sj ; rrfa£e_brishta is consistent with 
the lTikho nov fc Gal azutdinova ( 2005) analysis to larger 
scale heights. This suggests that their fit was dominated 
by the old stars, and not a younger population. However, 
the lack of an inflection in the surface density distribu- 
tion does not unambiguously rule out the presence of 
multiple components. 

The ubiquity of thick disks in galaxies has pre viously 
been proposed by iDalcanton &: Bernstein! 1)20021) based 
on color-gradients in edge-on disk galaxies. Our obser- 
vations confirm that the color-gradients (at least at the 



low-mass end) are the result of true differences in stellar 
populations. However, whether these gradients have an 
analagous formation mechanism to the Milky Way thick 
disk is not clear. 

One set of observations that reaches considerably 
de eper than thes e groun d-based observations is presented 
bv lZibetti et alJ l)2004(K who used stacked Sloan images 
to show that halos are common in late-type, edge-on 
galaxies. Their composite galaxy has a significantly 
wider field of view and poorer resolution than our ob- 
servations. They show that the best-fitting model to 
their data is a disk+halo model, with the disk component 
dominating out to roughly 10 exponential scale heights 
(^lOz^)- Their limited resolution and combination of 
a heterogenous sample of galaxies would likely prevent 
them from seeing the RGB components we see in our 
galaxies. However, the possible detection of the more ex- 
tended RGB components detected in IC 505 2, NGC 4144 
and NGC 5023 may be halos similar to the lZibetti et alJ 
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(l20f)l halo. 

4.3. Disk Heating 

The increase in zq seen in each galaxy between the 
MS, AGB and RGB populations (Fig. 7, Table 4) could 
result from a number of mechanisms, including vertical 
heating of a thin disk. Such a model would naturally 
produce the observed trend of older stellar populations 
having larger scale heights. In Figure 8 we plot the in- 
crease in scale height with mean stellar age for four fields 
that span the observed behaviour in our sample. We plot 
scale heights for the RGB and AGB, normalized by the 
MS scale heights, with the height of the symbols indi- 
cating the la uncertainties on zq. Note that when inter- 
preting our data in the context of disk heating models 
we are therefore implicitly assuming that the RGB and 
AGB stars were originally formed in a layer with a scale 
height comparable to that of the present day main se- 
quence stars. We assign characteristic ages to the RGB 
and AGB using the MATCH synthetic CMD tests for a 
constant SFR (§4.1.1). However, because the galaxies' 
actual star formation histories may differ significantly 
from the constant SFR assumed in Figure 4, we cannot 
assign a single age to each stellar population. Instead, we 
use the resulting age distributions to identify the 25th, 
50th (Median) and 75th percentile ages. The resulting 
age ranges are shown by the width of the individual boxes 
in Figure 8. However, note that the actual age of the 
population may lie entirely outside of the boxes, for ex- 
ample, if the RGB stars were all formed in a single burst 
12 Gyr ago. Thus, when interpeting Figure 8, one has 
substantial allowance in assigning an age. 

Overplotted on Figure 8 are dashed lines showing a 
range of power-law increases in the disk scale height z$ 
with time (zq oc t^ 13 ). For an isothermal sech 2 profile, 
the zo values are related to the vertical velocity dispersion 
{**): 

where E is the sur face density of the disk (Eq. 17 in 
Ivan der KrurdlT988|) . Studies of disk heating tradition- 
ally use power laws in the velocity dispersion, a z oc t~ a , 
and thus a — [3/2. Figure 8 therefore demonstrates that 
the vertical velocity dispersion of our galaxies has in- 
creased no faster than a = 0.15. More specifically, there 
are no characteristic ages that can be assigned to the 
AGB and RGB stars that yield heating rates greater 
than a = 0.15 (with the possible exception of the more 
massive, interacting galaxy NGC 4631), and thus this 
conclusion is robust even in light of our substantial age 
uncertainties. 

In contrast, the disk heating that has been observed in 
the Milky Way is comparatively rapid. The age-velocity 
dispersion relation (AVR) for Milky Way disk stars sug- 
gests that the vertical velocity dispersion increases with 
time with va lues of a ranging be tween 0.3 and 0. 6 (e.g. 
lWieledll"977t iBinnev et aLll2000t iNordstrom et aUl2004L 
see summary in Table 1 of Hanninen & Flynn 2002). In 
contrast, our limit of a < 0.15 is significantly smaller 
than the Milky Way value. These data immediately sug- 
gest that any disk heating in our low mass galaxies has 
been far less effective than in the Milky Way. Moreover, 
if some fraction of the extraplanar RGB stars are not due 



to disk heating, and are instead due to direct accretion 
or in situ formation at large scale heights, or if the RGB 
stars are weighted towards old ages (as we argue below 
in §5), then the actual rate of disk heating is even lower 
than suggested by Figure 8. 

There are several reasons why disk heating is ex- 
pected to be low for our sample galaxies. Within the 
Milky Way, the increase in vertical velocity dispersion 
with t ime is thought to be due t o scattering by spiral 
arms llBarbanis fe Woltieil |1967t iSellwood fc Carlberd 
Il984t ICarlberg fc Sellwoodl 119851) . by m olecular clouds 
(iSoitzer fc Schwarzschildlll951l). or both llCarlberglll987t 
[jenkjns^^Bninev[THj9(I I.TenkinsI Il992t IShapiroetal] 
120031 see also the review by Lacey 1991). However, our 
galaxies have sufficiently low masses and surface densi- 
ties that they are unlikely to be globally gray itationally 
unstable ijDalcanton et alJl2004t IVerde et al.ll2002|) and 
thus would not host strong spiral arms. Given that scat- 
tering by spiral arms seems to be the dominant heat- 
ing mechanism in the Milky Way (e.g. most recently 
iDe Simone et al JI200 4") . the absence of spiral arms alone 
should cause a drastic drop in heating rate down to 
a ~ 0.2 — 0.25, the expected value for heating by gian t 
molecular clouds alone (e.g. iHanninen" fc Flvnnl 12002). 
Likewise, the absence of strong dust lanes in these sys- 
tems and the results of >l4.4l both indicate that the cold 
molecular ISM is in a thicker layer than in the Milky 
Way. This large scale height for the cold ISM, and 
the general lack of molecular gas in low mass galaxies 
l|Young fc ScovillellT99ltfLerov et al.ll2005l) should there- 
fore further suppress the efficiency of disk heating. Fi- 
nally, the young stellar disks in our sample are apparently 
much thicker than in the Milky Way, wh ich could reduc e 
the efficiency of any heating mechanism l)Freemanlll991j) . 
iShapiro et al. ( 2003]) also argue for reduced disk heating 
in late-type galaxies based on the ratio of vertical to ra- 
dial velocity dispersions. However, their rationale for the 
observed trend is opposite from what we conclude from 
our data. 

As an aside, the low observed heating rates may 
provide strong constraints on cosmologically impor- 
tant sources o f disk heating incl uding late-time satel- 
lite accretion doiiinn et alJ I1993J) . massive black holes 
llLacev fc Ostrikerl 119851) or halo su bstructure (e.g. 
IHanninen fc Flvnd 120021 iBenson et alJl2004|) . However, 
the expected heating rates for such models have been 
calibrated for massive spiral disks, not the thicker, lower 
surface density galaxies studied here. 

4.4. Modelling Dust Effects on the Stellar Density 
Profiles 

Before continuing to explore the origin of extraplanar 
stars, we briefly examine the vertical distribution of the 
dust layer. At first glance, interpretation of the stel- 
lar density profiles near the midplane in Figure 5 might 
be somewhat confusing. If the dips in surface density 
are due to dust, why does the dust appear to affect the 
RGB and AGB stars more than the MS stars? We sug- 
gest this may occur because the dust layer is opaque 
near the midplane and is distributed with a scale height 
greater than or equal to the MS population, but less than 
the AGB/RGB populations. The MS stars we are see- 
ing would then lie entirely in front of an obscuring dust 
screen, while the AGB/RGB populations would have a 
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Fig. 8. — This figure shows the increase of scale height with time 
based on the stellar populations analysis described in §4.1-2. The 
AGB and RGB scale heights are shown for four galaxies spanning 
the range of zo . Each scale height is shown with a box, color coded 
by galaxy. The box width represents the 25th, 50th, and 75th 
percentile age range of the AGB or RGB population (see Fig. 4) and 
the height represents the ler errors in the determined zo (Fig. 7). 
The zq values were normalized by the zo of the MS population. 
The dashed lines show power-laws of the form zo cc t~@ . The RGB 
(upper) bar and AGB (lower) bar for any galaxy can be used to 
roughly constrain the disk heating rates. Because the precise star 
formation history of these components is not known, the width of 
the boxes is derived from a constant star formation rate model and 
thus represents the range of ages that contribute to the AGB and 
RGB population. 



significant population at large disk heights above where 
the galaxy becomes optically thin. The dip in their num- 
bers near the midplane is then explained because the op- 
tically thick dust layer obscures some fraction of the stars 
along the line of sight. 

To test this explanation, we built a simple 'toy model' 
galaxy with MS, RGB and AGB populations distributed 
as sech 2 profiles with the z$ as shown in Table 4. All com- 
ponents were given identical radial distributions with the 
K s band exponential scale length. The dust component 
was assumed to also follow a sech 2 profile with a variable 
scale height, zo.dust, and a radial distribution identical 
to that of the stars. For simplicity we assumed that the 
dust has no effect at an optical depth less than one, but 
is completely obscuring at greater optical depths. Thus, 
along a line of sight, the dust is completely transparent to 
t = 1, and completely opaque beyond. For each vertical 
position we integrated the dust component along the linc- 
of-sight until an optical depth of one was reached, which 
set the depth of the dust screen at that height. Stellar 
density profiles for the three separate populations were 
then created by totaling the number of stars in front of 
the dust screen at each height. We then normalized the 
stellar density profiles as in Figure 5. 

Figure 9 shows the resulting model of NGC 4144 for 
three values of zo,dust presented for comparison to the 
observations shown in Figure 5. In each case the amount 
of dust in the midplane is the same. The underlying 
values for z$ were adopted from Table 4 (374, 699, and 
934 pc for the MS, AGB and RGB respectively). The 
left panel shows the results for a dust layer whose scale 
height is narrower than all three stellar components. For 
this case, there is a pronounced dip in the surface density 
profile of all three components. In the middle panel, the 
value of zn.dust is between the zo values for the MS and 
the AGB/RGB populations. In this case there is a dip 



only in the AGB/RGB, because the height of the MS 
layer is entirely confined within the opaque dust layer, 
allowing only the unobscured stars on the near side of 
the galaxy to be detected. The right panel has a dust 
layer larger than the both the MS and AGB value and 
therefore a dip is seen only in the RGB component. The 
middle panel does a good job of qualitatively matching 
the observations for NGC 4144 in Figure 5. 

Referring back to Figure 5, we can see that the main 
sequence profile lacks a dip near the midplane for most 
of the fields, while the RGB has a dip in all cases within 
2-3 zi/2- This suggests that the dust in the galaxy is 
opaque below 2-3 21/2 and that it is distributed in a 
layer with thickness greater than or equal to the MS 
stars. Because the MS stars are already distributed in 
a thicker dis tribution than in t he M ilky Way, this result 
supports the iDalcanton et aD ((2004) finding that galax- 
ies with circular velocities below 120 km s _1 have large 
dust scale heights and do not form thin dust lanes. All of 
the galaxies presented here except NGC 4631 are below 
this circular velocity limit (Table 1). 

This model also suggests that although the depth to 
which we see in each galaxy is different at differing scale 
heights, it is the same for all the stellar populations at a 
single scale height. This validates the comparisons made 
in Fig. 6 between the ages of the stellar populations at 
different scale heights. 

Although this model matches the gross characteris- 
tics of many of the profiles shown in Fig. 5, it fails to 
fully explain their details. Most notably, in NGC 55 
and NGC 4244, the MS profiles are significantly lower 
than the best-fitting sech 2 function, contradicting the 
idea that all the stars we are seeing lie in front of a 
screen. This is likely the result of our model's lack of 
sophistication. Physically, it doesn't take into account 
the possibility of flares or changing scale lengths as a 
function of stellar population. In addition, it treats dust 
extinction in a very simplistic fashion. A more sophisti- 
ca ted treatment of the dust (such as the one presented 
bv lMatthews & Woodl l200l) is beyond the scope of this 
paper. The conclusions reached in this section should be 
considered tentative and will be tested in a later paper, 
in which we will present a more realistic dust model. 

5. METALLICITY DISTRIBUTION FUNCTIONS 
As has been shown previously (e.g. 



As has been show n previously (e.g. 
Da Costa fc Armandrol 119901 lArmandroff et all 119931 



Fravn fc G ilmorc 2002) , the color of the red giant branch 



near its tip can be used to constrain the metallicity of 
old stellar populations. Although reddening due to dust 
will prevent an accurate measurement of the metallicity 
of stars near the midplane, we can determine a rough 
metallicity for stars above the midplane where the effect 
of dust and the contamination from AGB stars are small 
(as shown in Figure 6, middle panel). 

Figure 10 shows a composite CMD of all high-latitude 
stars (above Az 1/2 ) in IC 5052, NGC 55, NGC 55-DISK, 
NGC 4144, NGC 4244, and NGC 5023. This disk height 
limit was chosen (1) to be well above any of the dips as- 
sociated with dust in Figure 5, (2) to dramatically reduce 
contamination by AGB and HeB stars that might inter- 
fere in the metallicity determination, and (3) to be where 
thick disk stars dominate in the Milky Way l)Chen et alJ 
12001(1 . Overlayed on the CMD are 10 Gyr old RGB 
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Fig. 10. — Composite CMD showing the red giant branches 
of IC 5052, NGC 55, NGC 55-DISK, NGC 4144, NGC 4244, 
and NGC 5023 at disk heights >4z\/2- The absolute magnitude 
Mf814W was determined using the TRGB distance modulus in Ta- 
ble 1. Solid lines show 10 Gyr old Padova models for the RGB 
with [Fe/H] (from left to right) of -2.3, -1.7, -1.3, -0.7, -0.4, and 
0.0. The peak of the stellar distributions fall between [Fe/H] of -1.3 
and -0.7. Dashed lines indicate the three bins used to determine 
the metallicity distribution functions shown in Fig. 11. 



isochrones ijGirardi et al J 121)0(11) at metallicities ranging 
from [Fe/H] = —2.3 to 0.0, with higher metallicities be- 
ing redder. Examination of Figure 10 shows that the 
peak of the distribution falls between the [Fe/H] = — 1.3 
and -0.7 lines. Roughly 13% of the stars fall bluewards of 
the [Fe/H] = —2.3 isochrone, probably due to a combina- 
tion of photometric error and the presence of AGB stars. 
Very few are redder than the [Fe/H] = —0.4 isochrone. 
Overall, Figure 10 indicates that most of the stars above 
4:Z 1 / 2 are moderately metal-poor. 

Improving on the metallicity determinations in Pa- 
per I, here we will derive metallicity distribution func- 
tions (MDFs) usi ng untransformed magnitudes and 
Padova isochrones l)Oirardi et alJl2f)0fy . To determine 
metallicity distribution functions for individual galaxies 
we binned the stars in up to three independent 0.2 mag- 
nitude wide bins centered on Mpsi4w of -3.5, -3.3, and - 



3.1. These are shown with dashed lines in Figure 10. The 
brightest bin was chosen to include stars as metal-rich as 
[Fe/H] = —0.4, which do not get brighter than Mp8i4w 
of -3.7. We considered only the magnitude bins that were 
above the 20% completeness cutoff for crowded regions 
(thus excluding the fainter bins in NGC 4144, see Fig- 
ure 1). This cut eliminated NGC 4631 and NGC 4631- 
DISK from the analysis, because none of the bins fell 
above their 20% completeness limits. We corrected the 
colors and magnitudes of each star for foreground redden- 
ing, but made no correction for completeness, since the 
change in completeness across the color-range in ques- 
tion is similar to the error that would be introduced by 
that correction. We then determined the metallicity of 
each star by linearly interpolating between the 10 Gyr 
isochrones in each bin. Since some of the stars blue- 
wards of the [Fe/H] = —2.3 isochrone may be RGB stars 
scattered to bluer colors by photometric error, we at- 
tempted to include these stars in the MDFs. In each 
bin, stars bluer than the [Fe/H] = —2.3 isochrone were 
given positive color shifts by multiplying a Gaussian ran- 
dom number by their error. This correction moved >50% 
of these stars redwards of the [Fe/H] = —2.3 isochrone, 
and resulted in a slight increase (~0.01 in the normal- 
ized units) in the MDF between a [Fe/H] of -1.3 and -2.3. 
Table 5 gives the resulting number of stars used in MDF 
determination, the median F606W-F814W error, and the 
peak and mean metallicities of the MDFs. 

Figure 11 shows the resulting MDFs for each field. The 
shaded regions show the total range in the MDFs as de- 
rived in the different magnitude bins, and indicate that 
our results are consistent among the different magnitude 
ranges. In general, the MDFs peak at metallicities [Fe/H] 
of -0.7 to -1.1 as expected from Figure 10, and have a tail 
of stars to low metallicities. 

Before analyzing the MDFs further, we note that 
there are several uncertainties in the detailed shapes of 
the metallicity distribution functions. First, they are 
b ased on isochrone models and not empirical data (as 
in lSaraiedini fc Van Duvnell200lD . Second, the color er- 
rors on the metal-poor end (Table 5) translate into a 
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TABLE 5 
Metallicity Determinations 



Field 


Paper I 


# of stars 


(7 


MDF Peak 


Mean 




[Fe/H] 




(F606W-F814W) 


[Fe/H] 


[Fe/H] 


IC 5052 


-1.22 


3770 


0.21 


-0.7 


-0.9 


NGC 55 


-1.41 


1068 


0.09 


-1.0 


-1.0 


NGC 55-DISK 


-1.62 


983 


0.06 


-1.1 


-1.2 


NGC 4144 


-1.50 


1773 


0.19 


-0.9 


-1.0 


NGC 4244 


-1.45 


1970 


0.11 


-0.9 


-1.0 


NGC 4631 


-1.25 










NGC 4631-DISK 


-1.54 










NGC 5023 


-1.71 


1119 


0.22 


-0.9 


-1.0 



Note. — The columns from left to right are: Field name, the metallicity estimated 
in Paper I, the number of stars used in the metallicity determination in this paper, 
the typical color error of those stars, the peak of the metallicity distribution function 
(from Fig. 11), and the mean metallicity. 



large error in metallicity. We estimated the errors as 
a function of metallicity by inserting stars of a specific 
metallicity/color, giving them appropriate color errors, 
and then determining the spread in the resulting metal- 
licity distribution. This procedure gave errors of 0.5-0.8 
dex at [Fe/H] = —2.3. However, the shape of the distri- 
butions are much more believable on the metal-rich end 
where the isochrones are well separated, giving errors of 
less than 0.1 dex at [Fe/H] = -0.4. At the peak of the 
MDF ([Fe/H] ~ -0.9), typical errors are 0.2 dex, sug- 
gesting that the peak metallicities derived here are fairly 
reliable. 

Comparing the peak metallicity of the extraplanar 
stars to known Milky Way populations, Figure 11 in- 
dicates that the metallicities of the extraplanar stars are 
a factor of ten times too high to be an alogs of the Milky 
Way' s stellar halo ([Fe/H] - -1.7, IWvse fc Gilrnorel 
1995). This result would not change even if all the stars 
bluewards of the [Fe/H] = — 2.3 isochrone are low metal- 
licity RGB stars. We note however, that the low metal- 
licity of the Milky Way halo may not be typical. The halo 
of Andromeda has been found to be much more metal- 
rich than in the Milky Way, with a pea k [Fe/H] ~ -0.6 
(Holl and et"al]ll996t Brown et al.ll2003j) . although there 
is difficulty ascribing these outer stars to a halo popu- 
lation per se, given M 31's complicated outer structure 
{Ferguson et alil2002D . 

Of all the Milky Way components, we find that the 
peak metallicities are most consistent with those of the 
metallicity of the Milky Way thi ck disk, which has [Fe/H ] 
- -0.8 based on F/G dwarfs JWvse fc Gilrnorel 1199^ . 
The extraplanar stars studied here are somewhat more 
metal poor than the Milky Way's thick disk (by up to 0.3 
dex). However, this offset may not be surprising given 
the lower mass of our galaxy sample (V c ~ 80 km s _1 vs 
V c ~ 220 km s" 1 for th e Milky Way). 

As in the Milky Way llWvse fc Gilmorell 995t lH avwood 
120011) . the metallicities of the extraplanar stars appear 
to be more metal poor than the thin, young, main se- 
quence population. Although dust prevents us from re- 
liably measuring the metallicity of stars near the mid- 
plane, we can estimate their metallicity using the current 
gas phase metallicity. NGC 55, the only galaxy in our 
sample with a gas phase abundance measure ment, has 
12 + log (O/H) = 8.32 at one disk scale length ijGarnettJ 



2002). This metallicity corresponds to [Fe/H] 0.6 

(assuming [Fe/O] = 0), which is 0.5 dex more metal rich 
than the extraplanar stars at a comparable radius. Other 
late- type disks with similar rotational velocities from the 
iGarnettl ((2002) compilation have comparable gas-phase 
metallicities, suggesting that the offset in metallicity be- 
tween the midplane and the extraplanar populations is 
likely to be systematic. 

Although they do not explicitly examine stars as a 
function of scale height, studies of metallicities in other 
galaxies using methods similar to ours also find broad 
agreement with the presence of an extended [Fe/H] ~ — 1 
population of RGB stars. The LMC (which has a mass 
similar to the galaxies in our sample) has a peak metal- 
licity distnbjrtionofJFe/H] ~ —0.6 for RGB stars in the 
disk l|Cole et alJ 120001) . Recent papers on the outer re- 
gions of M33 also find pe ak [Fe/H] values of -1.0 l)Davidgd 
200 31 iTiede et al.ll2004|) . Furthermore, a recent paper 
bv iDavidgg l)2005l) . derives a [Fe/H] of roughly -1 for 
NGC 55 using near-IR photometry of resolved extrapla- 
nar stars, closely matching our peak metallicity in Fig- 
ure 11. Our data and others therefore suggest the perva- 
sive presence of a significant [Fe/H] ~ — 1 old population 
in late- type galaxies. If our association of this popula- 
tion with a thicker disk is generally true in other galaxies, 
then it presents an attractive solution to the "G-dwarf" 
problem seen in the Milky Way by providing the neces- 
sary prompt initial enrich ment for stars in the thin disk 
()Truran fc CamerorJll97lT) . 

Overplotted on Figure 11 as dashed lines are the ex- 
pected metallicity distributions for c losed-box "s imple" 
chemical evolution models (Eq. 20 of lPagell|l997() scaled 
to the peaks of the MDFs. While the basic shape of these 
models are similar to our MDFs, there appears to be a 
deficit in some galaxies of stars at both low and high 
metallicities. A deficit at high metallicities is expected if 
star formation truncates before all the gas is consumed. 
Within the context of thick disk formation models, this 
truncation may occur if some of the gas reservoir that 
forms the extraplanar stars instead settles into the thin 
disk, if the extraplanar stars were heated from a previ- 
ously thin but gas-rich disk, or if the extraplanar stars 
were directly accreted from merging satellites that suf- 
fered from tidal stripping or supernova blowout. The ap- 
parent deficit of stars at low metallicities may be another 
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manifestation of the widespread G-dwarf problem. Thus, 
while the existence of a substantial population of stars 
at [Fe/H] ~ — 1 may help to solve the G-dwarf problem 
in the thin disk, it may have simply pushed the problem 
into a new component. The solution to the extraplanar 
G-dwarf problem will likely lie among the suite of pop- 
ular model s previously explored for the thin disk (see 
lPage]|ir99^ . However, some of the deficit of stars at low 
metallicities may also result from the photometric errors 
and methods used to construct the MDFs, as discussed 
above. 

Finally, we note that the peaks of the metallicity distri- 
bution functions given in Table 5 are significantly more 
metal-rich than our previous determination presented in 
Paper I. Rather than using native F606W— F814W col- 
ors, these earlier measurements a pplied the metallicity- 
color relation of iLee et alJ (|1093f) to the mean color of 
the giant branch transformed to the Johnson-Cousins fil- 
ter system. These previous values, reproduced in Ta- 
ble 5, range from [Fe/H] of -1.2 to -1.7, versus -0.7 to 
-1.1 in the present work. We believe that in addition 
to the magnitude transformation, the difference in de- 
rived metallicity results in part from the difference in 
binning (in [Fe/H] vs. F606W-F814W), such that the 
mean color does not correspond to the peak metallicity. 
This offset can be seen in the 0.1-0.2 dex offset between 
the mean and peak metallicity in Tabic 5. Furthermore, 
the Paper I determination also included stars at some- 
what lower disk heights, thereby increasing the number 
of AGB contaminants. We believe that the MDFs and 
their peak metallicity that we present here are more re- 
liable than the estimate given in Paper I. 

5.1. Vertical Metallicity & Color Gradients 

Models for the origin of extraplanar stars (i.e. disk 
heating, direct accretion, etc.) predict different degrees 
of variation in the stellar metallicity with height above 
the plane. To investigate the vertical variation of metal- 
licity with disk height z we have examined the median 
color of the RGB stars as a function of the height above 
the midplane. Figure 12 shows the median color of RGB 
stars between Mp8i4\v of -3.2 and -3.6 and redwards of 
F606W-F814W=0.7, binned by the scale height of the 
galaxies. Data are plotted where errors in the median 
color are <0.05 magnitudes. The hatched region at low 
disk heights shows where the effects of internal redden- 
ing may impact the colors of the stars. For the three 
fields with profiles extending beyond 2-3 kpc we bin the 
RGB stars at large disk heights in a single bin to reach 
adequate signal-to-noise in our measurement of the me- 
dian color, plotted as diamonds in Figure 12. We note 
that these are the same stars that comprise the possible 
extended components discussed in §4.2. 

Figure 12 demonstrates three main points. First, the 
color gradients in the galaxies are relatively small, indi- 
cating that the stars have nearly uniform metallicity with 
increasing distance above the plane, particularly at scale 
heights above the region potentially affected by dust. 
However, we note that the stars at very large radii (shown 
with diamonds) do tend towards bluer colors, possibly in- 
dicating the presence of a more metal-poor population at 
z > 10zi/2 (2-3 kpc). Second, the color gradients show 
no systematic trends, and are equally likely to be rising 
or falling. Finally, all the galaxies have very similar RGB 



colors (as demonstrated already in Figure 11). 

Our metallicity gradient results are consistent with 
previous ob servations of these and other low-mass spi- 
ral ga l axies (lDavidgeH2005l iTikhonov et al.ll2005l IMouldl 
12005ft . IMouldl (|2005r l used HST archival data to study 
the vertical properties of disks in four low-luminosity 
(My ~ —16) edge-on galaxies comparable to those stud- 
ied here. Using AGB, RGB, and red supergiant stars over 
a large range of magnitudes (—8.5 < Mi < —1.5), he cal- 
culates the mean colors up to 2 kpc from the plane. His 
main results are that there are slight or no color gradients 
as a function of disk height and that the metallicities of 
the stars at disk heights between 400 and 2000 pc are be- 
tween -0.8 and -1.0 in all four gal axies, in excellent agree- 
ment with what we find here. ITikhonov et al. (2005) 
also state that any metallicity gradient in NGC 55 is 
very small, in agreement with our observed lack of color- 
gr adients in Figure 12 . Recent simulations of thick disks 
bv lBrook et alJ l)2005f > also show a lack of any metallicity 
gradient with disk height. 

The lack of strong metallicit y gradi e nts ca n be ex- 
plained in a number of ways. IMouldl l|2005T) suggests 
the lack of metallicity gradients rules out dissipative and 
simple accretion models for thick disk formation, and 
favors a model in which thick disks form during interac- 
tions. We note that the lack of metallicity gradient may 
also be enhanced by an "age-bias" in the metallicity mea- 
surements. The mass of stars on the RGB changes from 
~2 M at an age of 1 Gyr to <1 M Q at 10 Gyr. Given 
the steepness of the IMF in this mass range, the RGB 
age distribution will be weighted towards older ages, and 
therefore to a more uniform metallicity. However, based 
on the constant star formation rate models in §4.1, it 
still appears that the RGB stars could span a wide range 
of ages, even in the presence of the expected age-bias. 
Therefore, the simplest interpretation of the lack of gra- 
dients in our data is that many of the RGB stars at disk 
heights above 42^/2 (~1 kpc) formed at a similar time 
and thus have comparable enrichment histories, elimi- 
nating any metallicity gradient. This scenario could be 
explained either by a sudden heating of the disk by in- 
teraction, or by accretion of gas-rich satellites which re- 
sulted in the formation of a thick component. N-body 
simulations have shown that early merging and accretion 
events can produce thick disks with old ages ijAbadi et alJ 
l2003HBrook et alJl2004 iBrook et al.ll2005[) . 

6. DISCUSSION & CONCLUSIONS 

The work presented here has identified a number of 
main observational results: 

• In low mass, late-type galaxies the thickness of a 
stellar population increases systematically with the 
age of the stars being studied. This behavior has 
been seen not just in all six of the galaxies studied 
here, but in all other HST studies of edge-on late- 
type galaxies, and in the Milky Way as well. The 
larger scale heights of older stellar populations is 
therefore likely to be a generic property of galaxy 
disks. 

• All of the studied galaxies show a clear interme- 
diate age (1 — 5 Gyr old) population whose scale 
height is intermediate between that of the young 
main sequence stars and the older red giants. 
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Fig. 11. — Metallicity distribution functions for stars at disk heights >4z 1 / 2 - Gray regions indicate the range of values obtained in 
different absolute magnitude bins. The bins used for this determination are 0.2 magnitudes wide and are shown in Fig. 10. Only bins above 
the magnitude limits in Table 3 were used. Dashed lines indicate closed-box chemical evolution models scaled to the peak of the MDF. 



• The metallicity of the dominant old stellar popu- 
lation has [Fe/H] ~ — 1, but shows little or no gra- 
dient between 3z 1 / 2 and 10zi/ 2 above the plane. 
Above this height (> 2 — 3 kpc), there are tenta- 
tive indications of decreasing metallicities, which 
may be associated with slight overdensities in the 
RGB surface density at similar distances above the 
midplane. 

• In the low mass galaxies studied here (V c ~ 70 — 
130 km sec" 1 ), the young stellar population is sys- 
tematically thicker than in the MW, and has a ver- 
tical scale height comparable to the thickness of the 
dust layer. This suggests that the cold ISM has a 
larger scale height in low mass galaxies, consistent 
with the lack of dust lanes observed in such systems 
Ipalcanton et al.ll2004|) . 

• The young and intermediate-age stellar popula- 
tions dominate the integrated near-infrared light 
of late-type low mass galaxies. 



We now interpret these observational facts in the con- 
text of disk formation models. First, taken at face value, 
the old RGB component's ~3:1 axial ratio and [Fe/H] 
~ — 1 metallicity suggest a close correspondence with the 
Milky Way's thick disk. However, each of our galaxies' 
scale heights steadily increase from the young main se- 
quence to the intermediate age AGB and the older RGB. 
The uniformity of this trend strongly suggests that our 
disks are not simply the superposition of two components 
(i.e. a thick and thin disk). Instead, the data require 
a more complex model incorporating some disk heating 
to explain the systematically larger scale height of the 
intermediate-age population. The necessary disk heating 
would also have affected any older population, and thus 
must make some contribution to the thicker population of 
RGB stars. The required amount of disk heating is much 
smaller than is seen in the Milky Way (§4.3), and could 
likely be provided by molecular clouds or minor mergers. 
The latter scenario is slightly favored by the large vari- 
ations in the apparent change of disk scale height with 
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Fig. 12. — The median RGB color as a function of disk height in 
each galaxy. This includes stars with F606W-F814W redder than 
0.7 and Mpgi4w between -3.2 and -3.6. The error on all plotted 
data is <0.05 magnitudes. The diamonds indicate the median color 
of all stars at large disk heights and is plotted only for the three 
galaxies with significant numbers of stars at large scale heights. 
The hatched area indicates the heights at which internal redden- 
ing may affect the stellar colors. The right ax is shows the mean 
RGB color of Padua isochrones iGirardi 2004) at three different 
metallicities. 



time (Figures 6 & 8) . Heating through satellite accretion 
and interaction could naturally produce these stochastic 
variations. However, numerical simulations of heating in 
diffuse low mass disks would be required to definitively 
constrain any of the above scenarios and to assess how 
significant a contribution heating may have made to the 
thicker RGB component. 

While the above argument strongly suggests that disk 
heating must play some role in the production of the 
extraplanar stars, the lack of a metallicity gradient in 
RGB stars at moderate disk heights (500 — 2000 pc, as 
shown here and in Mould 2005) suggests that steady 
disk heating cannot entirely explain the thickest com- 
ponent of old RGB stars. If the past star formation rate 
has been constant, then there is a significant overlap be- 
tween the stellar ages of the AGB and RGB regions of 
the color-magnitude diagram (Figure 4). A significant 
fraction of the RGB stars should therefore have smaller 
scale heights, younger ages, and enriched metallicities, 
and would thus produce a steady increase in metallic- 
ity towards the midplane for all but the most contrived 
scenarios. The most attractive explanation for the lack 
of metallicity gradient is that instead of a constant star 
formation rate, the majority of RGB stars at all scale 
heights must have formed early and with a well-mixed 
metallicity distribution. Such a population would dwarf 
any subsequent population of enriched RGB stars with 
lower scale heights. While steady dynamical heating 
could push this ancient population to larger scale heights, 
it could not simultaneously account for recent dynami- 
cal observations of counter-rotating disks at these disk 



heights in comparable galaxies l)Yoachim fc Dalcantonl 
2005). Taken together, these observations are better ex- 
plained by scenarios involving the fo rmation of a thick 
disk of stars in merger events (as in lAbadi et all feOQS: 
iBrook et all 120041 IBrook et al.ll2005l) . Overall, our re- 
sults require that some disk heating occurs at interme- 
diate ages (to puff up the AGB component), but that 
events at earlier times (interactions or mergers) created 
a majority of the RGB stars over a short timescale. 

Finally, we present tenuous evidence for an extended 
old component seen only at disk heights > 2-3 kpc. At 
large scale heights we see marginal overdensities of stars 
in the RGB profiles of Fig. 7. There also seems to be a 
reduction in the metallicity of R GB stars at this height 
(Fig. 12). In one of our galaxies. iTikhonov et alJ (|2005f) 
finds strong evidence for a very extended component of 
RGB stars extending from ^2-7 kpc. While this compo- 
nent appears to have an exponential distribution, its zq 
value of ^2 kpc (compared to a radial scale length of ~1 
kpc) strongly suggests it is not a disk. These extended 
components are detected at about the same height where 
th e halo becomes prom inent in the stacked Sloan images 
of iZibetti etafl ((2004). However, based on our observa- 
tions, we are unable to assess the properties or frequency 
of these components. 

The present-day structure of galaxy disks results from 
a complex mixing of effects and a full explanation re- 
quires detailed knowledge of the the star formation his- 
tory, merging events, and disk heating. Studies like the 
one presented here help to disentangle these effects and 
determine their relative importance as function of galaxy 
type and mass. This study also shows the promise that 
HST observations of resolved stars have for enabling the 
detailed analysis of low surface brightness stellar compo- 
nents in galaxies outside the Local Group. A comparison 
of our data with N-body simulations of low mass disk 
galaxies would assist in constraining disk heating and 
merging scenarios. Unfortunately, current simulations 
of disk galaxy formation have focused on massive galax- 
ies like the Milky Way l|Abadi et al J 12003 IBrook et all 
120041 IBrook et al. 2005). Also, deeper observations that 
fully resolve the red and blue horizontal branches would 
greatly improve constraints on the star formation histo- 
ries of these galaxies and improve our understanding of 
their structure. 
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